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Two programs have been written for the recently developed single- 
chip digital signal processor (dsp) integrated circuit that enable it to 
function as a tone generator in testing transmission systems. One 
program is based on the table look-up method and the other on the 
Maclaurin expansion method. A dsp tone generator based on the 
look-up method can generate up to 12 components and is suitable for 
all transmission testing applications. A generator based on the Ma- 
claurin expansion method is limited to less than four components 
and is particularly applicable in two-tone testing. 

I. INTRODUCTION 

A tone generator is required for a number of transmission system 
tests. The tones required include TOUCH-TONE® signaling, multi- 
frequency (mf) signaling, a milliwatt source, centralized automatic 
reporting on trunks (carot) responses and carot test tones. In addi- 
tion, a tone consisting of 21 components each having a settable phase 
and level is required for a fast Fourier transform (fft) -based system. 
In all cases, the tone is to be transmitted on a 4-kHz digital channel, 
with a sample every 125 ^is (Nyquist sampling rate of 8 kHz). Programs 
have been written for the dsp that enable it to function as a tone 
generator. The following two methods are used: (i) table look-up, and 
(ii) Maclaurin expansion. 

The table look-up method consists of storing in read only memory 
(rom) the trigometric values of sin(n6), where < n < N, and N is 
determined by both frequency granularity requirements and harmonic 
distortion considerations. At each sampling instant, the value of the 
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sample is taken from the appropriate location of the ROM table and 
scaled for its desired level. If more than one tone is desired, each 
component is independently determined and then all component val- 
ues added together to form the sample value. 
The Maclaurin expansion, 

cos(x) = l-|j+^--| [ +---, 

can be used to determine successive sample values by direct calcula- 
tion. The number of terms to be considered is a function of desired 
harmonic purity and the time allowed for the computation. 

II. TABLE LOOK-UP METHOD 

A program has been written for the dsp that uses the table look-up 
method to implement a tone generator. It produces up to six inde- 
pendently specified components in the 0- to 4-kHz range with an 8-bit, 
ji-255 encoded output. A diagnostic is included that checks most of the 
dsp features used by this routine. The program accepts successive 16- 
bit serial input instructions that: 

(i) Specify tone generation or a diagnostic. 

(ii) Specify the number of components, from 1 to 6, that comprise 
the tone. When the number of components specified is 0, quiet tone is 
generated. When 7 components are specified, milliwatt tone is gener- 
ated. 

(Hi) Specify the frequency of a component, from to 4 kHz, in 1/8- 
Hz steps. 

(iv) Specify the phase of a component from to 360° in 1/8° steps. 
(u) Specify the level of a component as a fraction of full output, 
from to 1 in steps as fine as 2" 15 . 

(vi) Steps 3, 4, and 5 are repeated for each component comprising 
the tone. 

Tone generation starts with a reset, except when phase continuity 
is required. The output may be changed to a new frequency without 
loss of samples or phase continuity when the information about the 
new frequency is presented to the input buffer. 

The diagnostic ends with the SO-bit set to 1, if all tests pass. 



2. 1 Algorithm 

A full 2tt sine table, with 512 16-bit entries is used. A tt/2 sine table 
could be used but would require keeping track of quadrants and signs, 
which takes time and reduces the number of tones that could be 
generated. 

The table entries are so arranged that the entry for 0° is at an 
address called &TABLE. Two hundred fifty-five successive locations 
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above this contain the positive half of the sine wave trignometric 
entries. Two hundred and fifty-six entries below this address contain 
the negative going entries. Two additional entries, corresponding to 
the first nonzero positive entry, at location (&TABLE - 257) and at 
(&TABLE + 256) are also included. Hence, the rom table is as shown 
in Fig. 1. 

Table accesses above and below the normal range may occur because 
of the effect of rounding. Electing truncation in the dsp in the initial- 
ization of this process would avoid the need for adding these memory 
locations, but their presence assures continuity of sample values. 

The frequency of a component is determined by the number of table 
entries A<£ stepped per sample, if the samples being generated at a 
particular frequency are 

y(nT) = A sin(27r/hT), 

where f is the desired frequency, T is the sampling interval, and n 
takes on integer values, then the phase increment between successive 
values is 

2irfT = 2tt j, 

fa 



where f B ■ — = 8000 Hz is the sampling frequency. 
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Fig. 1— Read only memory table. 
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Since a full 2-n sine wave table is present in ROM and is represented 
by 512 entries, the phase difference between successive entries is 
2tt/512. The number of entries, A<#>, to be stepped per sample is, 
therefore, 

2.//8000 
v 277/512 

To minimize processing time, modulo 512 arithmetic is used so that 
no sign or table limit checking is required. Incrementing around 
the table is automatic. To accomplish this, the table is entered at 
&TABLE + N+ where N+ is modulo 512 and is stored in register y as 
shown in Fig. 2. As desired, N+ is never larger than 255 and its sign 
alternates. Assuming an initial phase of 0, 2V + starts at zero and is then 
incremented by A$. Thus, for p passes 

ty- E tt*. 

Given f, A<f> is obtained in modulo 512 by transferring as shown in Fig. 
3. 

To generate a frequency close to the desired value, it is important to 
have as much precision as possible in A<f>. In this case, after w register 
truncation, it is 11 bits. The precision also has a marked influence on 
harmonic distortion.* The harmonic distortion products are shown to 
be 50 dB below the fundamental for this arrangement. 

2.2 Phase continuity 

To change the frequency, A<f> must be changed, but phase continuity 
can be preserved. In one sample interval, the new frequency is fetched 
and the new A<f>, subsequently calculated, is placed into the appropriate 
random access memory (ram) location. The calculations then proceed 
from the previously accumulated phase. 

2.3 Flow chart 

A flow chart of the dsp program is given in Fig. 4. In essence, the 
overall flow breaks down into the following areas shown in the chart: 

(la) Input a data control word, the frequencies, phases and levels, or 

(lb) Input a diagnostic control word and run. The output is a pass 
or fail indication. 

Because of the dynamic nature of the dsp ram, it is necessary during 
the input data routine, to refresh the ram locations. 

Once the input data process is completed, one of three paths is 
possible: 

(2a) Calculate sample values for 1 «£ n < 6 tones. 



* W. N. Fabricius, unpublished work. 
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Fig. 2 — Modulo 512 stored in y register. 



(2b) Calculate samples for quiet tone, n = 0. 
(2c) Calculate samples for a dBm, milliwatt, n = 7. 
(3a) Continue calculating with phase continuity. 
(3b) Continue with old input data. 

2.4 Experimental results 

The spectrum of the dsp output signal after decoding is given in 
Figs. 5, 6, and 7. Calculations, based on the work of Fabricius, indicate 
that the harmonic production is mainly because of the jU-255 format. 

Some synthesized waveforms are represented by the waveforms in 
Figs. 8 and 9. 

The measured phase jitter is 0.65° at 1000 Hz and 0.75° at 1004 Hz. 
The am jitter experienced is less than 0.1 dB. 

The linearity of the output signal fell well within the ^-255 format 
requirements. (For this experiment, a linear output was delivered by 
the program.) 

2.5 Extensions 

The program, without jii-255 encoding, has generated 13 components 
within a 125-/us sampling interval. A total of 21 components could be 
generated by synchronizing two dsps, one producing 10 and the other 
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II. The first dsp would serially transmit its 20-bit sample value, after 
generation, to the second dsp. This would take 8 /is at a 2.5-Mb/s rate, 
and could be accomplished during the time the second dsp was 
computing its 11th component. The second dsp would then take this 
value in its input buffer and add it to its sample value, convert the 
result to /x-255, and transmit it. 

III. MACLAURIN EXPANSION 

A program has been written for the dsp that can evaluate a truncated 
Maclaurin expansion in less than 30 /is. Frequency resolution is not 
limited by quantizing as in the table look-up method, but by the 
accumulator and the product registers in the dsp. 

3. 1 Algorithm 

Each component is generated by a truncated Maclaurin series ex- 
pansion. The choice of the number of terms included in the approxi- 
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1 KILOHERTZ 



981 HERTZ AT -61 dBmO 




1 KILOHERTZ 
Fig. 5— Tone spectrum— low level (981 Hertz). 



1 662 THE BELL SYSTEM TECHNICAL JOURNAL, SEPTEMBER 1 981 



1962 HERTZ AT dBm 



1962 HERTZ AT -60 dBm 




1 KILOHERTZ 



1 KILOHERTZ 
1962 HERTZ AT -68.9 dBm 




1 KILOHERTZ 

Fig. 6— Tone spectrum at 1962 Hertz. 

mation and the point about which the series is expanded all affect the 
harmonic distortion. 

Consider the following Maclaurin series expansions around zero, 

q3 fl5 nl 

sin(0) = 6- — + + ... 

3! 5! 7! 

/l2 /)4 /)6 

The accuracy of the approximation, of course, becomes better when 
more and more terms of the expansion are used. The harmonic 
distortion because of series truncation of the cosine series is examined 
in the Appendix. Let P{6) represent the portion of the series that is 
retained, and R(6) represent the remaining terms so that 

cos d = P(d) +R(d), 

where both of these functions are considered defined only in the 
interval 

-ir/2 < 6 < it/2. 
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1 KILOHERTZ 

Fig. 7— Tone spectrum— high level (981 Hertz). 

They are replicated, with sign reversals, in successive intervals of -n 

forming a wave periodic in the interval 2-n. An upper bound on the 

harmonic distortion, because of truncation after the nth term of a 

Maclaurin expansion, is 

T (2n )! 
Distortion < = 20 log 2 , 

(TT/Z) 
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Fig. 8 — Time domain waveforms. 
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Fig. 9 — Spectrum of time waveforms. 

Table I lists this upper bound for values of n. This is a mathematical 
upper bound and ignores the effects of quantizing distortion. 

The harmonics generated due to the /i-255 format generally are in 
the —42 dB range. The total signal to distortion is in the 36-dB range. 
Thus, n = A should be used for the Maclaurin expansion since it 
produces distortion components less than those generated by the /x- 
255 format. 

By a similar analysis, it can be shown that the sine series should 
also be terminated after four terms so that the harmonic distortion 
upper bound is —78.9 dB. The Maclaurin cosine series is chosen instead 

Table I — Harmonic 
distortion 



Upper Bound 



-33.61 dB 
-60.73 dB 
-91.97 dB 
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of the sine series, since its highest term contains fl 6 instead of 6 1 and, 
thus, less programming is required. 

The program starts at — tt/2 and increments in steps of <£ = 277/T, 
where fis the desired frequency and T is the interval between samples. 
The computed function should remain positive until after n samples, 
when the accrued phase 6„ is such that 9„ = —it/2 + rufy and exceeds 
7r/2. At this time, tt is subtracted from the accrued phase yielding a 
new value 6 n such that —tt/2 <d n <0 and the process starts again. The 
computed value is now made negative. The negative sign prefixes each 
sample until, again, the accrued phase exceeds +77/2. Once again, it is 
subtracted from the accrued phase. The sign of each computed point 
becomes positive, and the algorithm cycle is complete. 

An expansion around tt/4 would result in less error in the resulting 
approximation at the end points, but it would require more program- 
ming to fix the sign. The errors for an expansion around zero are less 
than those attributable to the ]U-255 operation and, thus, it is not 
important to reduce them. 

3.2 Flow diagram 

A simplified diagram, Fig. 10, covers the generation of one compo- 
nent. Note that the loop counter (lc) serves as a flip-flop for deter- 
mining the sign of the result. Suppose lc = and tt/2 has been 
exceeded. The right-hand branch decision (lc^O) point in Fig. 8 causes 
lc to be decremented to —1. During the next pass when the left-hand 
decision (lc^O) is reached, the sign of the result is changed. This is 
maintained until the accumulated phase would again exceed +77/2, 
when traversing the right-hand decision point forces lc = 0. Now, the 
computed result is positive and the process repeats. 

3.3 Experimental results 

Figure 11 shows the spectrum of the dsp output after decoding. The 
harmonics are approximately 42 dB below the fundamental as pre- 
dicted for a jU-255 decoder. Note, however, the differences in the 
spectrums for frequencies of 500 and 502 Hz. Since 500 Hz is a rational 
submultiple of 8 kHz, all samples per sine wave cycle are repeated in 
subsequent cycles. There is essentially no quantizing noise, i.e., the 
sample values are completely periodic. The pronounced spectral lines 
demonstrate this purity. For the 502-Hz case, there is no periodicity in 
the samples for successive cycles and, hence, quantizing noise exists. 
The peak value of the spectral harmonics are reduced and the noise 
floor is raised. 

IV. CONCLUSIONS 

A dsp tone generator using the table look-up method is able to 
generate up to 12 components (or 21 components for 2-coupled gen- 
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Fig. 10 — Maclaurin expansion flow diagram. 

erators) assuming a 4-kHz, /x-255 channel with an output sample every 
125 /is. Harmonic distortion is substantially below that inherent in the 
channel. Such a versatile generator is suitable for all transmission 
system testing applications. 

Given a suitable decoder, a dsp tone generator can also be used to 
create any analog wave representable by 13 spectral lines restricted to 
the 4-kHz band. A sine wave component can be generated every 17.75 
jits and a 26-kHz sine wave can be generated. By suitable program 
changes, it seems possible to produce a 50-kHz sine wave. 

A dsp tone generator using the Maclaurin expansion method is 
limited to less than four components. It requires less ROM than the 
table look-up method and, hence, is applicable where such a program 
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Fig. 11 — Maclaurin spectrums. 

would be co-resident with other measurement routines in a single dsp. 
Where 2-tone testing is required, e.g., envelope delay distortion, slope- 
sag intermodulation distortion, etc., such a generator is especially 
applicable. 
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APPENDIX 

This analysis is directed to specifying an upper bound on the 
harmonic content of a truncated Maclaurin expansion. 
Given the Maclaurin expansion 

02 04 06 08 02 ( „-l» 

C0S , = 1 __ + ___ + __ + _____.., (1) 

consider it as being made up of a finite polynomial P(6) consisting of 
the first n terms of the above expansion and the remaining portion of 
the series R(6), i.e., 

cos 6 « P(6) + R(9). 

Consider now that a periodic function be generated such that 



and 



P N (6)=P(6) for ~\^Q<\ 



Pn(0) = elsewhere. 



(2) 
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Similarly, 



R N (6)=R(6) for -|s0<|, 



and 

R N (6) = elsewhere. (2a) 

By a replicating procedure, then 

cos0= £ (-l) n P N (0 - nir) + J (-1)"R,,(0 - nir). (3) 
The portion of eq. 3, 

(X 

cos s £ (-1) b Pjv(0 - nir), 

It— 00 

is the approximation used for sample generation in the Maclaurin 
series program given in this paper. 

Because of the definitions in eqs. (2) and (2a), each sum specified in 
eq. 3 is periodic with period 2-n and, hence, can be expanded into a 
Fourier series. Let pk and r k be the Fourier series coefficients defined 

by 

i (-l) n P N (0 - nir) = l p# M , (4) 

It 00 k— 00 

l (-l) n R N (6 - nir) = i ne ike (5) 

It— 00 k— 00 



and 



where 



and 



cos0= E pke*°+ £ ne' kB , (6) 

A— —00 A— -co 



p k = JL I V (-1)"P„(0 - n7r)e-' M dB 

27T 



^ = 7rl £ (-1)"«n(0 - mr)e- iM dd. 
2tt J_ ff „-_. 



(8) 



Interchanging the order of integration and summation in eq. 7, then 
Pk = ^ £ (-1)" I Pn(0 - nir)e-*°d6. (9) 

277 „ — oo I _ 
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Changing the variable of integration y = 6 — nir modifies this equation 
to be 



ITT n =- x 



-■nin-U 



P N (y)e- iky dy, 



—ir(n + l) 



and further simplifies to 



1 



^77 n — „ 



1 P N (y)e- iky dy. 

-w(n+l) 



(10) 



(11) 



Because of the definition of Pbiiy), the integral only exists for n = 1, 
-1, 0; 



Pk = 



r/2 



P N {y)e- iky dy + (-l) (1+k) P N (y) e - iky dy 

LJ-t/2 J-rr/2 



n/2 



+ (-1)- (1+A) P N {y)e- iky dy 



Combining these terms results in 



p k = — [1 + (-l)« + »] PN(^)e-'^. 

2,7 J- w /2 



r/2 



. (12) 



(13) 



When k is even, pk = 0, therefore, only odd harmonics exist in the 
periodic function given by eq. 4. Since Pn(Q) is an even function, it 
follows that 



Pa = - 



»/2 



P N (6) cos kddO k odd. 



(14) 



By a completely similar development for r b , and substituting into 
eq. 6, 



cos 6 = 2 



£ p* cos &0 + £ r * cos £0 

f/odd kodu 



From this equation, it is apparent that 

Pk = — r* and k ¥* 1 



(15) 



(16) 



since only the fundamental exists on the left-hand side of the equation. 

Given the properties of R(9) such that each term is positive and 

smaller than the preceding term, then one can state that the magnitude 

of the first term of R(0) is greater than \R(6) |. Hence, if P{9) contains 
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n terms then 



riln 

h.-\ Rm (17) 



The stipulation that the terms are successively smaller implies for the 
m th and m lh + 1 term that 

. „ , , 2m(2m - 1) 
1 < [^ 2 ""- 1, /(2(m - l))!]/[^ 2 ' n /(2m!)] = ^ . 

Since this must hold for all -it/2 < < tt/2, then 

2m(2m - 1) > I ^ 

and, hence, for this inequality to hold 

m > 1 applies. 

Equation 17, therefore, applies for maximum 0, after the second term, 
and certainly then 

O/2) 2 7(2n)!]>|#(0)| for n > 2. (18) 

Parsevals theorem states that 



2 

' ' J (-D'lMf-iMr) 



277 



d0= E Hb. d9) 



/.'odd 



Since i?(0), and hence, #n(0), are bounded, as has been stated in eq. 
(18), then certainly 

[(7r/2) 2 V(2n)!] 2 > 2 HI- J pl (20) 

*=3 A-3 

Aodd Aodd 

Therefore, the first term ofR(0), evaluated at jt/2, represents an upper 
bound on any harmonic component of the replicated Pn(0) i.e., the 
truncated Maclaurin expansion. The harmonics of 

(Replicated P N (0)) < 20 log [(2n)!/(w/2) 2 "]. (21) 
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